Cosmological consequences of generalised RS II braneworlds 
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We discuss certain features of cosmology in a generalised RS II braneworld scenario. In this 
scenario, the bulk is given by a Schwarzschild-anti de Sitter or a Vaidya-anti de Sitter black hole 
in which an FRW brane is consistently embedded, resulting in modifications of the 4-dimensional 
Friedmann equations. We analyse how the scenario can be visualised and discuss the significance of 
each term in these modified equations both for early time and for late time cosmology. We further 
analyse the perturbation equations, based on Newtonian as well as relativistic perturbations and 
show that the scenario has the potentiality to explain structure formation by the "Weyl fluid" arising 
from embedding geometry. The results thus obtained are confronted with observations as well. 
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INTRODUCTION 

The Randall-Sundrum braneworld model fl| provides 
us with a unique mechanism of solving the hierarchy 
problem and obtaining an effective 4-dimensional New- 
ton's law of gravity valid at large as well as small length 
scales, by the method of embedding our 4D universe in a 
higher dimensional anti-de Sitter (AdSs) space. Of late, 
a more general, and geometrically rich, picture of obtain- 
ing an effective 4D theory of gravity, by projecting the 
bulk field equations onto the brane, has been formulated 
by Shiromizu et al [2]. In this scenario, the bulk space- 
time is not necessarily AdSs, rather a generalised version 
of it. The so-called effective Einstein equation on the 
brane, thus obtained, not only gives important insight 
to the bulk-brane interplay but also raises the possibil- 
ity of explaining the gravitational phenomena of the 4D 
universe from a broader perspective. Since the formal- 
ism was brought forth, it has been applied to explain 
several gravitational aspects, leading to interesting re- 
sults which justify further investigation on the prospects 
of braneworlds. A thorough discussion of these issues is 
available in the review of S. Kar 3] in this volume. 

The cosmological aspects of this theory were developed 
at a rather later stage, due to the complicated equations 
arising from the extra terms involved in the theory, which 
makes the theory difficult to probe at a first go. The ba- 
sic aim of this article is to provide a brief review of what 
has been done in the cosmological context of braneworlds 
so far and show that, in spite of great complications in- 
volved, the cosmological consequences of this scenario re- 
sults in interesting physics which is worth investigating 
further. 



FRW BRANE, SCH-ADS BULK 

In the braneworld scenario, the bulk may be either 
empty with only a bulk cosmological constant, or it may 
consist of non-standard model fields minimally or non- 
minimally coupled to gravity or to brane matter. For 
empty bulk, the bulk metric in which an FRW brane can 
be consistently embedded, is given by a 5-dimensional 
Schwarzschild-Anti de Sitter (Sch-AdSs) black hole [i,[i| 
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where S3 is the 3-space and the function f{r) is given by 

(2) 
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where /c(= 0, ±1) is the spatial curvature. Here m is 
a constant which is identical to the 5D analogue of the 
Schwarzschild mass of the bulk black hole. 

An important point to note here is that the bulk met- 
ric presented here is asymptotically anti-de Sitter. When 
the mass of the 5D black hole vanishes, this metric can 
be recast in the warped RS II form by coordinate trans- 
formation. Therefore, the scenario is a generalisation of 
the RS II braneworlds, which reflects cosmology on the 
brane. In this sense, this setup may be called generalised 
RS II braneworlds. 

The Sch-AdSs black hole provides a novel way of visu- 
alising cosmological phenomena on the 4D universe. In 
this scenario, the brane is moving in the bulk, with its ra- 
dial trajectory given by r{T), where t is the proper time 
on the brane. With u^ = {i, r) the velocity vector on 
the brane, the normalisation condition is given by 
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g^i.u'^u'' = ~-fP 



(3) 



and the unit normal vector n^ , satisfying the orthonor- 
mahty conditions, is given by 



Further, the brane matter conservation equation 
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With this, we can get the induced metric on the brane 

ds^ = -dT^ +r^{T)dnl (5) 

which is, clearly, an FRW metric on the brane, with the 
scale factor a{T) being identified with the radial trajec- 
tory r{T). So, the interpretation here is as follows: What 
we call an expanding universe is visualised by a brane- 
based observer only. For a bulk-based observer located 
somewhere outside the bulk black hole horizon, this ex- 
pansion is exactly identical to the movement of the brane 
along the radial trajectory of the black hole. 

Given this situation, the important question is: How 
do the Friedmann equations on the brane look like? The 
review of S. Kar in this volume 3] discusses in details that 
the 4D Einstein equation on the brane is a generalisation 
of the standard Einstein equation. The so-called Effective 
Einstein Equation is given by (4I 

(^fii' = ~J^9iii' "f" >^i^ i-iv + K^Of^i, — t^i, (b) 

where the 4D cosmological constant and coupling con- 
stant are related to their 5D counterparts by 
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5p^ is the quadratic contribution from brane energy- 
momentum tensor 



^'^ ^ 12 



^r^aT;+ir„^r"V.-^rV (8) 

and Sfii, is the bulk Weyl tensor, projected onto the 
brane, which is given by 

(9) 



Thus, along with the usual contributions from the cos- 
mological constant and matter on the 4D universe, the 
EEE contains two additional terms : 

• Sp,^ : the local correction from brane matter 

• £fj_i, : the nonlocal correction from bulk geometry. 

For a perfect fluid on the brane, using the symmetry 
properties, these braneworld corrections can be added 
up to the usual brane energy-momentum tensor. The 
effective energy density, pressure, momentum density and 
anisotropic stress thus obtained are, respectively, given 
byi: 
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leads to the conservation equation for the Weyl term (via 
4D Bianchi identity on the brane VG^jy — 0) 



4V = 
a 



Hence the Weyl density p* behaves as 



(15) 
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where C is a constant which is basically the rescaled bulk 
black hole mass. 

Further, for a bulk compatible to FRW geometry on 
the brane, q* = = ■Tr*^,. With the use of the effec- 
tive Einstein equation ^ and the effective quantities 
obtained in equations (fTO]) and pT|) for a perfect fluid 
on the brane, the Friedmann equation and the covariant 
Raychaudhuri equation turn out to be 
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Obviously, these equations are quite different from the 
standard 4D equations due to the presence of the 
quadratic terms and the Weyl terms. 

Let us now explore the roles of these additional terms 
on 4D cosmology. The quadratic term iS^;^ plays a signif- 
icant role in the early universe {p^ ^ A{,). For example, 
this leads to a faster Hubble expansion at high energies 
and a more strongly damped evolution of the inflaton 
field [6| . Thus the brane universe inflates at a much faster 
rate than what is expected from standard cosmology. 
Further, with this term, the modified equations make 
it possible to explain the inflationary scenario which is 
driven not by a inflaton field on the brane (i.e., not by 
any 4D field as required in standard cosmology) but by a 
dilaton field in the bulk [7|. However, the quadratic term 
is negligibly small at late times as p^ <C A{, > (lOOGeF)'*. 
Hence, its role is relevant so far as only the early universe 
is concerned. 

On the other hand, the role of E^v is to supply an addi- 
tional perfect fluidlike effect to the actual on-brane per- 
fect fluid (since, from the previous arguments, 9u = = 
7r*j,). The so-called Weyl fluid arises as the tidal effect 
of the 5D black hole in the bulk and its density contri- 
bution p* is related to the mass of the bulk black hole. 
Hence, in order to have a realistic contribution from the 
Weyl term on the brane, we need the black hole mass to 
be positive, so that p* > 0. Since for a vacuum bulk, Eq 
P^ reveals that the Weyl fluid is strictly radiation-like. 



it does not play any significant role in late time cosmolo- 
gies. It can, at best, slightly modify the standard per- 
turbative analysis. Its role has been extensively studied 
for metric-basedperturbations [8| , density perturbations 
on large scales [9|, curvature perturbations [10| and the 



Sachs- Wolfe effect 
perturbations [1^ 
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vector perturbations [12IJ . tensor 
[i|. In all the 



and CMB anisotropics 
cases, the effect has been found to be slightly enhanced 
from the standard analysis. However, we shall show in 
the subsequent discussions that the Weyl fluid can have 
a crucial role in late time cosmologies as well when the 
bulk is not necessarily empty. 



FRW BRANE, VAIDYA-ADS BULK 

Till now we have discussed brane cosmology with 
empty bulk. However, as already mentioned, the bulk 
may, in principle contain non-standard model fields, for 
which the above situation is further generalised. In this 
case, the bulk metric, for which the FRW geometry on the 
brane is recovered, is given by a radiative 5D Vaidya anti- 
de-Sitter black hole (VAdSs). In terms of transformed 
(null) coordinate v = t + J dr/ f, the bulk metric can be 
written as 



dSl^ 



~f{r, v) dv^ + 2dr dv + r^dT,l 



(19) 



where S3 is the 3-space. For a spatially flat brane, the 
function /(r, v) is given by 



fir, v) 



m{v) 
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(20) 



with the length scale I related to the bulk (negative) cos- 
mological constant by A5 = —6/P and m{v) is the vari- 
able mass of the Vaidya black hole. 

In the case of such a radiative bulk, the effective Ein- 
stein equation on the brane is further modified to 15|, Il6l | 



Clearly, in this scenario, the extra terms come from 
three sectors : 



S. 



liu 



a quadratic term from brane energy- 
momentum tensor 

• E^i, : a geometric term from the bulk Weyl tensor 
projected onto the brane 

• J-'fj^i, : a term involving the brane-projection of the 
bulk energy-density. 

The combined effect of the last two terms is related to 
the sumtotal of the mass of bulk black hole and the radi- 
ation field, and this is now the Weyl fluid that supplies 
an additional perfect fluid-like effect to the usual brane 
perfect fluid. 



Here the exclusive contribution from the bulk matter 
is given by the energy-momentum tensor of a null dust 



T^bulk 
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(22) 



where qa are now the ingoing null vectors and ip oc dm/ dv 
is the rate of ingoing radial flow to the bulk black hole, 
ra{v) being the resultant of the masses of VAdSs black 
hole and the radiation field. This type of bulk can ex- 
change energy with the brane along the radial direction. 
There are extensive study in the literature on the energy- 
exchange between bulk and brane. (see, for example, 
[iSl llq). Consequently, the brane matter conservation 
equation is modified to 



p + i-ip + p) 
a 



-24) 



(23) 



which now contains an effect of bulk matter onto the 
brane and guarantees the null flow from the radiative 
bulk black hole. Also, the Bianchi identity on the brane 
V'^G^jy — leads to the equation governing the evolution 
of the Weyl fluid 
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(24) 



The difference of the above equation with Eq p3|) for 
empty bulk is remarkable. For a bulk with matter, the 
coupling term involving ip determines the nature of the 
Weyl fluid. A general expression for this is now given by 

M 



p = 



C(r) 



(25) 



T being the proper time on the brane. Clearly, when the 
bulk is not strictly matter-free, then the Weyl parameter 
C{t) is no longer a constant, and consequently, p* no 
longer behaves like radiation. It is only when the bulk 
is empty, we get back the radiation-like behaviour of the 
Weyl fluid from Eq ([24]) . 

It is obvious from the preceding discussion that the 
Vaidya- AdS bulk scenario is so far the most generalised 
description of the braneworlds for which cosmology of 
the 4D world is relevant. Thus, this is the most general 
generalised RS II braneworld scenario. 

The components of the effective energy-momentum 
tensor T°|J are now given by [l7| 



„eff 



P 
P 



2Xb 

P 

2Xb 



C{t) 



{P + '^P) 



C{r) 
3a4 



(26) 

(27) 



As before, for a Vaidya- AdS'5 bulk compatible to FRW 
geometry on the brane, both q'^^ and ttJ;^ vanish. We 
are thus left with a perfect fluid like effect on the brane. 



constituted of brane and bulk matter-energy and bulk 
geometry, with an evolving Weyl fluid. 

The final modifications to the Friedmann equation and 
the covariant Raychaudhuri equation on the brane are 
due to the collective efi^ect of the embedding geometry, 
brane matter and bulk matter. In terms of the brane and 
bulk quantities, these generalised equations read [l8| 
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It is worth mentioning that these equations are rad- 
ically different from those for empty bulk due to the 
presence of an evolving Weyl parameter C{t). So, they 
should, in principle, give rise to new physics on the brane. 



"NEWTONIAN" PERTURBATIONS 

Given the modified Friedmann equations, we will now 
engage ourselves in studying cosmological perturbations 
on the brane. Following Newtonian analysis of perturba- 
tions from gravitational instability, we will demonstrate 
that the growth of perturbations of the Weyl fluid can 
take care of the fluctuations required by the inhomoge- 
neous local universe, and thus, have the potentiality to 
explain structure formation without the need for dark 
matter |17| . 

The equations of hydrodynamics now involve the 
quadratic brane correction and the Weyl fluid correction 
to the brane perfect fluid. In terms of the effective quan- 
tities, these equations are given by 
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where if^ is the velocity field in the effective perfect 
fluid. It should be noted that the term $°* is not the 
usual Newtonian potential but the effective gravitational 
potential which is the resultant effect of the brane and 
bulk parameters in the form of effective quantities. 

Let us now consider small perturbation about the ini- 
tial unperturbed effective quantities. Perturbation on the 
effective density (the so-called effective density contrast) 
is given by 



/«(ir, t)=p^V)(1+^'=''(^, r)) 



(33) 



whereas the perturbation in the effective gravitational 
potential is 



$<="(^, r) 
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(34) 



where p^^{t) and <I>q*'' are respectively the unperturbed 
effective density and effective potential and S°^ and (j)°^ 
are their corresponding fluctuations. For completion, we 
mention here that a perturbation in the Weyl fluid would 
mean a perturbation on the bulk geometry. There is a 
possibility that this may destabilise the brane itself. In 
this context, however, we assume that the brane remains 
stable even after such perturbations, which can only be 
guaranteed if one analyses the effects of perturbations on 
the full 5D bulk metric. 

Expressing in terms of comoving coordinates v = 
a r + u'^^ and neglecting terms of second or higher order 
in the equations we arrive at the simplifled perturbation 
equations 
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The above set of equations have a unique solution given 
by 



5^*^(ir, r)=V<5f(r)e»^-^ 



(38) 



From now on, we shall express the perturbations in terms 
of the inverse Fourier transform of the above: 



5f{r)^^J5^\-t, r)e-^-^d3 



(39) 



For a barotropic fluid, the effective pressure is a function 
of the effective density only. Hence, equations ([35| - (|37p . 
with the Fourier mode solution, transform into a linear 
second order differential equation for the perturbation 

5f 
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where c^g is the square of the effective sound speed which 
is given by 
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The main job now is to analyse the second order differ- 
ential equation (|40)l for the effective perturbation. The 
perturbations of the effective fluid will grow at late times 
and will account for the required amount of gravitational 
instability only if the Weyl density redshifts more slowly 
than baryonic matter density, so that it can dominate 
over matter at late times. Recall that for a Sch-AdS 



bulk (empty bulk), the Weyl fluid is strictly radiation- 
like, hence it redshifts at a faster rate than baryonic mat- 
ter, eventually being negligible at late times. But for the 
VAdS bulk (radiative bulk), the nature and evolution of 
the Weyl fluid is governed by Eq ([M]) . This equation can 
be written conveniently as 



a 



(42) 



where Q is a coupling term which is determined by the 
projected bulk energy density ijj. Since the bulk infor- 
mations are not known ab initio, this coupling term is 
given by a physically reasonable and consistent ansatz. 
Certain earlier attempt of choosing an ansatz are listed 
in [lal2li|24|- We note here that, in order to study the 
evolution of the Weyl fluid as function of the scale factor, 
the general ansatz for the coupling term should be of the 
form ^17,] 



Q = aHp* 



(43) 



with a > 0. For this type of ansatz, the Weyl fluid 
behaves as 



1 



p oc 



l(4- 



(44) 



Consequently, Eq ([^5]) reveals that the Weyl parameter 
is given by 



CiT)^Coa''iT) 



(45) 



where Cq is its initial value at the matter-dominated 
epoch. Note that the Weyl fluid is strictly radiation- 
like only if a = 0, for which we recover the Sch-AdS bulk 
scenario, which is a special case of this generalised de- 
scription. So, in general, the nature of the Weyl fluid 
depends on the coupling strength a. In order that the 
Weyl fluid dominates over matter the coupling strength 
a should have a value within the range 1 < a < 4, as ob- 
tained from this "Newtonian" analysis. Later, we shall 
put more stringent bounds on the value of this parame- 
ter from relativistic analysis as well as from observational 
ground. 

Let us now get back to the the perturbation equation 
(|^0| of the effective density which is a sum-total of three 
quantities. As discussed before, the quadratic term being 
negligible at late times, the effective density at late times 
is practically given by 



^(6) 



(46) 



where p^^' is the baryonic density. Thus, along with the 
usual matter density, here we have an additional (Weyl) 
density contributing to the total density that governs the 
perturbation equation. This Weyl density, being geomet- 
ric, is essentially non-baryonic. Consequently, we decom- 
pose Eq (|40p by separating the baryonic (matter) part 



from the non-baryonic (Weyl) part, which results in two 
separate equations, one each for each of the fluids 
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where S^''' and S* are the fluctuations of baryonic matter 
and Weyl fluid respectively. In the above equations, the 
term involving sound speed has been neglected as we are 
interested only in the growing mode fluctuations. With 
57^'') ^ Q*, the relevant growing mode solution for Eq 
can be expressed as a function of the redshift as 



S*{z)^6*{0)il + zy 



(49) 



which, when put back into the fluctuation equation (|47|) 
of baryonic density, gives 



dr^ 



2-4— =4^Gp*r(0)(l + z)-i (50) 

a dr 



Since the late time behaviour of the expansion of the 
universe in RS II is the same as the standard cosmological 
solution for the scale factor 15|, |20| , for a spatially flat 
(fc — 0) brane, we have the scale factor at late time 

2/3{w + l) 
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Considering ri*°* 
simplifles to 
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ri* w 1 at present time, Eq ((50|) 
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which readily gives a solution of the form 17 1 
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where the scale factor is related to the redshift function 
hy a (X {1 + z)^^. 

Eq (|53p reveals that at a redshift close to z^, the bary- 
onic fluctuation 5'^'^^ almost vanishes but the Weyl fluc- 
tuation 6* still remains finite. This implies that even 
if the baryonic fluctuation is very small at a redshift of 
ZN ~ 1000, as confirmed by CMB data [l9|, the fluctua- 
tions of the Weyl fluid had a flnite amplitude during that 
time. On the other hand, at a redshift much less than zn 
the baryonic matter fluctuations are of equal amplitude 
as the Weyl fluid fluctuations. This is precisely what is 
required to explain the formation of structures we see 
today. 

The bottomlinc here is that in this analysis, we get 
a result which resembles standard cosmological relation, 
but has an essential distinction that nowhere we have to 



introduce dark matter by hand. Rather, here the fluc- 
tuations are governed by the Weyl density fluctuation 
d* which is a product of the embedding geometry via a 
modified Einstein equation in the braneworld scenario. 

Another point is worth mentioning here. The effective 
equation of state parameter 

,,,-ff _ ^ ^ P^'^ + P^'Hp^'^ + 2p(''))/2A, + C(r)/3a4 
ff ~ 
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reduces to (in this context of matter-dominated era) 
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which bears significant difference from the equation of 
state of cold dark matter (u; = 0). Thus, the theory gives 
rise to a geometric fiuid which is very different from dark 
matter in origin and nature but has the potentiahty to 
play the role of dark matter in cosmological context. 



RELATIVISTIC PERTURBATIONS 

Let us now proceed further to study relativistic pertur- 
bations on the brane based on the above analysis. The 
scenario, in brief, is as follows: Here the cosmological dy- 
namics is governed by a two-fluid system, out of which 
one is the material fiuid p^ ' , which is the ordinary (bary- 
onic) matter, and another a geometric fluid p* , termed as 
the Weyl fluid. These two fiuids interact and exchange 
energy between them in such a way that the total (ef- 
fective) density on the brane is given by p = p^^' + p* . 

M 

Given this scenario, the comoving fractional gradients 
of density and expansion are expressed, following usual 
general relativity, as 
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Further, the conservation equations (P5)l and ([M)l can 
be written composedly as 



P« 



+ e(pW+p«) = /« 



(59) 



where a superscript (i) denotes the quantities for the i- 
th fluid and I^^' is the corresponding interaction term. 
Written explicitly in the bulk-brane scenario discussed 
here, the interaction terms are: 



jib) 



-2^j 



I* = 2ip- 



2 /Kgy 

3 V k/ 



a 



(60) 
(61) 



In what follows we shall restrict ourselves to the dis- 
cussion of the Einstein-de Sitter brane universe for which 
il* = 1, ilj^ = 0. Here, with the above notations, the lin- 
earised evolution equations turn out to be 
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We shall now try to express the above equations in 
terms of covariant quantities. These density perturba- 
tions are governed by the fluctuation of the following co- 
variant projections 

AW = a Z)^ A|;) , A = a L*^ A^, Z = a Dt'Zf, (64) 

The covariant density perturbation equations on the 
brane, when expressed in terms of the above covariant 
quantities, turn out to be 
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Obviously, the equations are too complicated and it is 
almost impossible to have a possible solution from these 
complicated equations a priori. However, the equations 
turn out to be tractable if we incorporate certain simpli- 
fications following physical arguments. We have seen in 
the case of Newtonian analysis that for a radiative bulk 
the Weyl fluid, in general, evolves as 



3* = Cna-^^-°''> 



(67) 



with the parameter a in the range 1 < a < 4. We now 
consider V' to be a function of time only. This essentially 
means that we are considering only the time-evolution 
of the bulk-brane scenario, which is relevant for pertur- 
bation analysis. We further assume that the energy ex- 
change between the two fluids is in equilibrium. This 
basically describes the late time behaviour. Once again 
this is consistent so far as density perturbations are con- 
cerned. Consequently, the Weyl fluid now behaves as 



p (X a 



-3/2 



(68) 



with the parameter a = |. This readily suggests that 
the Weyl fluid redshifts more slowly than ordinary mat- 
ter and hence, can dominate over matter at late times, 
reflecting one of the fundamental properties of dark mat- 
ter. This also provides a more stringent bound for the 
value of a from theoretical ground alone. Later, we shall 
confront this value with observations. 

With the above results, the evolution equation for the 
Weyl fluid at late times is radically simplifled by using 
(syn << A*, which can be recast in convenient form as 



A* 



^A- 
t 



A* =0 



(69) 



where the constants A^ B, C are readily determinable. 
The above equation for A* turns out to be somewhat 



tractable. One of its solutions is given by [23 1 



A* 



1 2 2 Bessel/ 



y/l-2A + A^+ 4C, 2VBVt 



(70) 

We know that the Bessel function is a growing function. 
Therefore, the evolution equation for the Weyl fluid, in- 
deed, shows a growing mode solution, which is required 
to explain the growth of perturbations at late times. 

Thus, even the relativistic perturbations show that so 
far as theoretical results are concerned, the Weyl fluid 
can explain structure formation. 



CONFRONTATION WITH OBSERVATIONS 

An important issue is to confront this theoretical model 
with observations. Recently there has been some study 
in this direction [2J], although an extensive study still 
remains as an open issue. In this section, we shall very 
briefly mention what has been done in this direction so 
far. 

In terms of dimensionless parameters 

„ _ A K^po ^ 2Co _ k'pI 



3Hf " 3Hi' 



4-arr2- 



6Ai?2 



with the total density satisfying the critical value 



E' 



"tot — / "i — 1 



(71) 



(72) 



the Friedmann equation (including the contribution from 
brane tension) can be re-written as 

The luminosity distance for FRW branes is now given by 

dL {Z) = ^i±^ C 



For r^A ^^ relevant for late-time cosmology, this integral 
can be evaluated as 



lAA* 



= d^CDM 



i i^ J J: 



(74) 



■'°- [OAa6-|-Opaga3+0,ag"°a°+2+nAa|;] 



where /, is a function having elliptic integrals of 1st and 
2nd kind, which have exact analytical expressions. 

Comparing this with standard ACDM scenario, one 
finds that a certain amount of Weyl fluid with 2 < a < 3 
is in nice agreement with SNe data 2j|. 

It is important to mention here that in the relativistic 
analysis we have found that a = |, which falls in this 
region, and hence, the braneworld model of perturbations 
is so far an observationally accepted model. However, as 
already mentioned, an extensive study in this direction 
remains as a very crucial open issue, which, we hope, will 
be addressed in details in near future. 



SUMMARY AND OUTLOOK 

We have discussed certain features of cosmology in a 
generalised RS II braneworld scenario, where the bulk 
is either a Schwarzschild-anti de Sitter or a radiative 
Vaidya-anti de Sitter black hole. We have shown that 
the theory leads to a modified version of the brane Fried- 
mann equations. Specifically, the local corrections to the 
Friedmann equations are manifest via a quadratic contri- 
bution from the brane perfect fluid whereas the nonlocal 
corrections supply a Weyl fluid which arises as an effect 
of the bulk-brane geometry. We have investigated for the 
role of each of the terms for early times as well as for late 
time cosmologies. Further, we have shown that the Weyl 
fluid plays a crucial role in late time cosmology, for the 
most general bulk metric. We have demonstrated that 
its presence radically changes the perturbation equations, 
and have discussed some of the implications of fluctua- 
tions involving it. The fluctuations are found to grow at 
late times and thus, may take care of the large amount of 
inhomogeneities observed in the local universe. We also 
mentioned some observable sides of this model. 

Some future directions related to branewolrd cosmol- 
ogy involve both theoretical and observational aspects. 
In the theoretical side, a thorough study of different pa- 
rameters related to cosmological perturbation need to be 
performed. In the observational side, confronting this 
model with observations in more details will test the the- 
ory in a more conclusive way. Last, but not the least, 
the expansion history of the universe from these modi- 
fled equations can be progressed further, which may lead 
to more interesting results and may even be investigated 
to find out any possible link with observations. 
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